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ABSTRACT 

I demonstrate that one can detect pixel gravitational microlensing events at 
the rate ~ 180 events yr~^ and that some fraction of events (~ 15) will be good 
enough to measure time scales at the 20% level if the experiment is carried 
out toward the central M31 bulge region with a dedicated 1 m telescope and a 
2K X 2K CCD. I also show that the event rate could be significantly increased 
(up to a factor 6.5) if more aggressive observation strategies are adopted. These 
aggressive strategies include multiple observation sites, using a bigger telescope, 
and reduction of seeing differences among images. A few years of observation 
could provide a data set big enough to extract important information about the 
M31 bulge such as the optical depth distribution, the mass spectrum of M31 
MACHOs, and the luminosity function, which currently cannot be obtained by 
other methods. Additionally, I find that the extinction of M31 bulge light by 
dust in the disk is not important and can be approximated by the model in 
which only 20% of light in V band is extinguished over most of the bulge area. 

Subject headings: astrometry - dark matter - gravitational lensing - M31 
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1. Introduction 



Pixel lensing, gravitational microlensing of unresolved stars, was first proposed by 
Crotts (1992) and Baillon et al. (1993). Classical lensing experiments are restricted to 
monitoring only fields of resolvable stars, e.g., toward the Galactic bulge and Magellanic 
Clouds (Alcock et al. 1995; Udalski et al. 1995; Ansari et al. 1995; Alard et al. 1995). By 
contrast, pixel lensing can overcome this restriction and thus can extend our ability to 
probe both Galactic and extragalactic Massive Compact Objects (MACHOs). Currently, 
both groups have demonstrated the technique by detecting various types of variable stars 
(Ansari et al. 1996; Tomaney & Crotts 1996) 

For classical lensing, the light curve is obtained by subtracting the reference flux, B, 
from the amplified fiux, F: 

Fi{A-l) = F-B- B = Fi + B,^„ (1.1) 

where Fj is the flux of the lensed star after (or before) being lensed and i?res = J2j^i is 
the residual flux from stars not participating in the event. The ampliflcation is related to 
the lensing parameters by 

= 4)1/2 ; ^' = ^\t-t,f+p\ (1.2) 

where co = is the inverse Einstein ring crossing time, /3 is the angular impact parameter 
in units of angular Einstein ring radius, 9^, and to is the time of maximum amplification. 
The angular Einstein ring radius is related to the physical Einstein ring radius, Ve, and to 
the physical parameters of the MACHO by 



/4GM DoiDi 



1/2 



(1.3) 



where Doi, D\s, and Dos are the distances between the observer, source, and lens, and M is 
the mass of the lens. On the other hand, it is impossible to measure directly the absolute 
values of F and B for the case of pixel lensing because stars cannot be resolved. However, 
the light curve, Fi{A — 1), still can be obtained by subtracting the reference image from 
the image in which a lensing event is occurring. 

Due to its characteristic differences from classical microlensing observations, additional 
analysis is required to estimate the pixel lensing event rate. For classical lensing, an event 
with modest maximum amplification (e.g., ^max ~ 2) can often be detected with high 
signal-to-noise ratio, S/N, when observed for a few minutes on a 1 m telescope. This is 
partly because source stars are nearby, and thus bright, and partly because the residual 
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flux, Bj-cs, is moderate or negligible. Note, however, that Bj-^s is not exactly zero due to 
the blending of star light in dense fields. For the case of pixel lensing, however, stars are 
located relatively far away, and thus yield low signal. Additionally, the signal suffers from 
significant noise from background stars that are not being lensed. Therefore, the event 
detection rate per observed star is significantly lower than that from classical lensing and 
thus events can be detected only when either the amplification or the intrinsic flux is 
very high. However, the loss of detection rate per star is compensated by the fact that 
many more stars are monitored. Therefore, one can obtain an equivalent number of total 
detectable events to that of classical lensing experiments. 

In this paper, I show that one can detect pixel lensing events at a rate 
r ~ 180 events yr~^ with a dedicated 1 m telescope and a 2K x 2K CCD. Among 
the events, ~ 15 events have high enough S/N to measure the time scale at the 20% level. 
The rate can be increased by up to a factor ~ 7 by employing a few improved observational 
strategies, and thus ~ 100 events with measurable time scale could be detected per year. 
These strategies are discussed in § 7. Gould (1996) developed a general formalism for the 
analysis of pixel lensing. The reader is referred to this work for a better understanding of 
the formalism that I use without detailed derivation. Comparison of the formalism with 
that developed by Colley (1995) is also recommended. 



2. Formalism 



The signal-to-noise ratio during a single observation is given by 

iS/Nf,=^^[A{U)-lf, (2.1) 

where F is the flux of the lensed star, S is the surface brightness of observed area, fipgf 
is the angular area of the point spread function (PSF) , tcyc is the observation duty cycle, 
a is the mean photon detection rate averaged over the duty cycle. If the observations 
are carried out for 4 hrs every night and the camera can detect e photons hr~^, then 
tcyc = 1 day and a = ex (iexp/24) photons hr""*^, where texp is the mean exposure time per 
day. Individual images obtained during the whole event are combined, yielding the total 
signal-to-noise ratio 

is/Nr = ^^t,yoY:[Mu)-i?. (2.2) 
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By converting the sum fcyc J^ii-^iU) — 1]^ into an integral /f^ dt[A{t) — 1]^, one can obtain 
a compact form of S/N, 

where the factor is evaluated and shown in Figure 1 of Gould (1996). 

Events are detected either when the star is bright or the amplification is high. This 
can be quantified by setting a maximum impact parameter /3max within which a lensing 
event with flux F is detected for a certain threshold signal-to-noise ratio {S/N),^^^^. The 
maximum impact parameter is then found by solving equation (2.3), i.e., /3max is the 
solution of the equation 



/?max V TTF^a 

where O = is the inverse of median Einstein time scale. Once /3max is determined, the 
event rate per angular area dfl is computed by 

r/r 2 

— =CnX]ro/3max(i^,)<^(i^,); Tq = -a>T, (2.5) 

j 

where r is the optical depth, (piFj) is the luminosity function (hereafter LF) of M31 bulge 
stars, and is a normalization factor, = 



3. Optical Depth Distribution 



For the computation of the event rate, it is necessary to compute the optical 
depth [see eq. (2.5)]. The computation of optical depth toward M31 requires some 
geometrical transformations since M31 is highly inclined, i ~ 77° (Kent 1989). I define 
the intrinsic coordinates, {x,y,z); r'^ = + y^, that are centered at the center of M31 
and r and z are measured along and normal to the disk plane, respectively. These 
coordinates are related to another set of coordinates, {x',y',z'), which have the same 
center but are measured along the projected major, minor axis, and the line-of-sight, by, 
f = [(y'cosi + z'sini)^ -|- x'^]^/^, z = \y' smi — z'cosz|. The adopted distance to M31 is 
c^M3i = 770 kpc. 

There arc several populations of lenses when stars in the M31 bulge are monitored. 
The first is MACHOs in the Galactic halo. The mean optical depth to a single M31 star 
lensed by the Galactic halo MACHOs is 

i-rrG /■'^Msi 

TMW,halo = — 9— / "-DolPMW.halo 

Jo 



iDoi)D, 



3.1 
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wherc D = Do\D\s/ Do^- For the Galactic halo matter distribution, I adopt a modified 
isothermal sphere with a core radius rc,Mw; 

(r\-l ^0(<MW + ^o)/«MW + r < 200 kpc , . 

10, r > 200 kpc, 

where the normalization is po = 7.9 x 10~^ M© pc~^, i?o = 8 kpc is the solar galactocentric 

distance, and the core radius is rc,MW = 2 kpc (Bahcall, Schmidt, & Soneira 1983). The 
value of the Galactic core radius is uncertain, but it hardly affects the optical depth 
determination because M31 is located ~ 119° away from the Galactic center. With this 
model, I find TMw.haio = 4.4 x 10""'''. 

The second lens population is the MACHOs in the M31 halo. Since the rotation curve 
of M31 is flat at large radii up to r = 30 kpc from the center (Kent 1989a), the M31 halo 
is modeled by a standard isothermal model with a core radius Tc m31 ; 



-c,M31 /„ o\ 



where the core radius is determined from the relation 



With the total core mass of Mc^msi = -^M3i,disk + Mivisi.buige = 13.2 x lO-*^*^ Mq 
(see below) and the M31 rotation velocity of Vc,M3i = 240 km s~^, I find that 
''c,M3i = (2/7r)(G-Mc^M3i/^^c M3i) ~ ^-^ ^P^- Then the optical depth contributed by M31 
halo lenses is computed similarly to equation (3.1), and results in rM3i,haio = 1.91 x 10~^. 

Finally, the rest of the optical depth toward the M31 bulge is contributed by the M31 
disk+bulge self-lensing events. This contribution of optical depth is computed by 



TM31 = — 2~ / dDosPM3l{Dos) / dDolPM3l{Dol)D / dDosPM3l{Dos) 



(3.5) 



where pM3i = PM3i,disk + PM3i,buige, PM3i,disk and PM3i,buige are the M31 disk and bulge 
densities, and di = (Imsi — 30 kpc and ^2 = c/m31 + 30 kpc are the lower and upper 
boundaries of mass distribution along the line of sight. The M31 bulge is modeled as an 
oblate spheroid with its (unnormalized) matter density as a function of semimajor axis 
provided by Kent (1989b). The best-fitting value of the axis ratio, c/a = 0.75, is found 
by fitting the computed surface brightness (corrected for / extinction) to a F-band image 
of M31 bulge (see § 5 for the observations). In the central region where the M31 bulge 
mass dominates, the rotation velocity is ~ 275 km s~^. Thus I normalize the total M31 
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bulge iiicLSs within t — 4 kpc by -^^m31, bulge — 

(275 km s V^''c,Mw)^MMW,bulgc, where 
I'c.MW = 220 km s^^ and Mmw, bulge = 2 x 10^*^ Mq are the adopted rotation speed of the 
Galaxy and Galactic bulge mass within r = 4 kpc (Zhao, Spergel, &: Rich 1995). This 
normalization procedure results in -MM3i,buige = 47r/o^ drr 2pM3i,buige = 4.9 x lO^^ Mq, 
which is in reasonable agreement with the determined value of 4.0 x 10^° M© by Kent 
(1989b). The M31 bulge is cut off at r ~ 8 kpc. The M31 disk is modeled by the double 
exponential disk, i.e., 

PMU,,M .) = ^ exp (--^) exp (- A) , (3.6) 

where Sq = 280 Mq pc"^ is the normalization, and the radial and vertical scale heights 
are hz = 400 pc and hji = 6.4 kpc, respectively (Gilmore et al. 1990; Gould 1994). This 
results in MM3i,disk(-R < 30 kpc) = 8.3 x 10^° Mq. 

The resulting optical depth distribution caused by M31 population (M31 disk+bulge 
self-lensing plus M31 halo) is shown as a contour map in Figure 1. The total optical depth 
both by M31 and Galactic halo lenses are obtained by adding TMw,halo = 4.4 x 10^'^ to 
the values marked on the contour. The contour levels are marked on the map in units of 
10~^. I note that the determined optical depth is subject to many uncertainties, especially 
in the mass distributions. For example, Braun (1991) has argued that there is no need to 
invoke a massive dark halo component. Instead of a halo, his model has higher disk and 
bulge densities, c.f., MM3i,buigc = 7.8 x 10^° Mq and MM3i,disk = 12.2 x 10^° Mq within 
30 kpc from the center. Pixel lensing experiments may resolve the conflicts in the matter 
density distribution in M31 from the comparison between the theoretical predictions and 
future observational result. 



4. Time-scale Distribution 



Not only the optical depth but also the median time scale is required for the 
computation of the event rate. The time scale is determined by the combination of both 
the transverse speed, v, and the Einstein ring radius of the lens, i.e., te = r^jv. The 
transverse velocity is related to the lens geometry and velocity of the source, Vg, the lens, 
vi, and observer, Vq by 

As ' 



V 



Vl 



Dr.. 



+ Vo 



A 



(4.1) 



For M31-M31 self-lensing events, one can approximate equation (4.1) as 

V = Vl — Vg since -Dqs ~ -DqI- Then the rms transverse velocity is found by 
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(■y2^i/2 _ ^^2 _|_ ^2 _ 2(^vi • Vg))"*^/^ = \/2(t's) because the cross-term (v; • v^) = 0. Here 
Vs ~ \/2o"bulgc) a-iid cJbulgc ~ 156 km s^^ is the adopted velocity dispersion of tire M31 
bulge (Lawrie 1983). The Einstein ring radius is related to the location of the lens system 
by Te oc D^/^, assuming all lenses have the same mass. In the computation, I assume 
M = 0.3 M0 because the most probable lens population is low mass stars. For other 
masses oc M^/^. The probability of finding an event located at a certain location D is 
found by 

P{D) cx / ' o?L'osP(^os) / °' dDMDoi)D'/\ (4.2) 
Then the time scale probability distribution is computed from P{D) by 



f{Q = / > 

Jo 



dt'J 



' , {4GMD)^/^ 

to 



e 

VC 



PiD). (4.3) 



I find te by taking the median value of f{te)- The resulting median time scale ie for M31 
self-lensing events as a function of position {x',y') shown in units of days in Figure 2 as a 
contour map. For other population events, i.e.. Galactic halo and M31 halo events, the 
time scale distributions are similar to that of M31 self-lensing events (Han & Gould 1996b). 



5. Extinction by the Dust in the M31 Disk 



The light from the M31 bulge suffers some extinction by dust in the disk because half 
the M31 bulge lies in front of the M31 disk while the other half lies behind it. Hence the 
bulge light from behind the disk is affected by the dust. By experiencing extinction, the 
apparent LF of the stars behind the disk will be different from that of stars in front, and 
thus r is affected by the extinction [see eq. (2.5)]. Since the disk is tilted, the total light of 
stars influenced by the extinction differs from position to position. 

To quantify the amount of extinction, CCD images in V and A bands were obtained on 
the night of 1995 August 4 using Tek2K detector on 1-m United Sates Naval Observatory 
Ritchey Chretien Aplanatic telescope, with a plate scale of 0."68 pixel"^. The A band filter 
is centered at 6900il and has a band width of AA = 2500A. Bias and sky subtraction of 
the images were carried out using the interactive data reduction package IRAF. 

I find that the disk of M31 is quite transparent. Let Sback ^ind Sfj-ont be the surface 
brightnesses contributed by stars located behind and in front of the disk, respectively. 
Then the flux ratio between two bands is given by 

n'{x',y') = l^''-^'^i''l^''-\ (5.1) 

■^^,front + JA^A,hack 
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where / is the fraction of light that passes through the disk in each band. Let 
7^^j,j^ = Sy/S/i be the (instrumental) ratio of flux in the two bands in the absence of 
extinction. This quantity can be determined by measuring the photon counts along the far 
minor axis where nearly all of the bulge is in front of the disk and where, consequently, 
extinction has almost no effect on the color. Then, one can rewrite equation (5.1) as 

^ ^ T^'{x\y') ^ Sfront + /ySback 2) 
^crit ^front + fA^ back 

The fraction is related to /y by 

/A = /f; q = AA/Av^Q.%, (5.3) 

which implies more absorption in V band, i.e., TZ < 1.0. The surface brightnesses, Eback 
and Sfrontj are computed from the mass model in § 3 and the color ratio Tl{x',y') is 
obtained from the images of the M31 bulge. Then fy is determined by solving equation 
(5.2) numerically using the relation (5.3), i.e., by finding the solution of the equation. 

In general, equation (5.4) has two solutions, e.g., for 7?. = 1 both /y = (totally opaque) 

and /y = 1 (totally transparent) solutions are possible. For the actual data, one set 
of solutions has (/y) ~ 0.1, while the other set has (/y) ~ 0.8. If /y = 0.1, the M31 
bulge surface brightness distribution would be extremely asymmetric with respect to the 
photometric major axis because most of the light would be extinguished on the near side 
while only a small fraction would be lost on the far side due to the high inclination of 
M31. In fact, the isophotes are close to symmetric. The mean extinction is therefore 
fv = 0.8. Because the disk is optically thin, the overall event rate depends only on 
the mean extinction. I therefore adopt the simplified model of uniform extinction with 
/y ~ 0.8, which corresponds to // ~ 0.88. 



6. Realistic Event Rate Estimate 

6.1. Additional Restrictions 

The event rate formula in equation (2.5) is derived under the assumption that the 
dominant form of noise comes from source stars, i.e., the Poisson photon noise N^. 
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However, except in the central region, ~ (6.5 x 6.5) arcmin^, most of the M31 bulge region 
has a surface brightness less than the sky brightness, which is (/) = 19.5 mag arcsec"^ 
averaged over the phases of the Moon. Therefore, the noise due to the sky-background, 
Ng-ky, becomes important. In addition, the differences in PSF from image to image taken at 
different times cause extra noise (Gould 1996), which I denote Npsi- I define the additional 
factors in noise due to sky and PSF variation by 



Vsky 



1/2 



Vpsf 



(6.1.1) 



The values of these factors are computed by 

1/2 



??sky 



/ ^M31 + ^sky \ 
V SM31 / 



Vpsi 



^mi ^exp f '^^see 



1/2 



(6.1.2) 



where Sgky and Smsi are the surface brightnesses of the sky and M31, e^j is the rate of 
photon detection at the fluctuation magnitude, m/ = 23.2 mag for M31 (Tonry, 1991), 
S9see/dsee is the fractional difference in seeing between images, and ndiv is the number 

of sub-exposures in each exposure time tcxp- By dividing the total exposure time into 
smaller exposures, the noise due to the time variation of PSF decreases because A'psf is 
averaged out, while remains the same. As a result of the additional noise, the total 
noise, Ntot = (-^^ + -^gly + -^psf)^^^' increased relative to the pure M31 photon noise by 
a factor, 

= (^sky + ^psf)'/'- (6.1.3) 



Ntot 



The net influence of additional noise on the event rate is that the required {S/N)Tnin 
increases by the factor r]. 

However, not all the events that satisfy the strict (5'/A^)min test qualify as detectable 
events. Additional tests await for them. One of these tests comes from the duration of 
an event: if the event lasts for very short time, the event could not be detected however 
high the S/N is, because the event continues less than the observation interval (1 day). 
This restriction of time scale is quantified in the computation by counting only events 
with tcff ^ 0.5 day, where the effective duration time is defined by t^s = Pmaxtc (Gould 
1995). Some of these short events might occur during the observation which continues 
for iexp = 4 hrs each night in average. However, they would not be detected because of 



another restriction discussed below. 



The other restriction comes from the finite size of the source star. If the source size is 
bigger than the Einstein ring, 6^ > 0^-, the point-source amplification formula in equation 
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(1.2) is no longer valid at the central region of the light curve, i.e., when 

where = {Do\/ Dos){R* / v) ~ R^/v is the crossing time over the star with a radius i?*- 
Then the maximum amplification will be overestimated by equation (1.2) and the event 
with 0.5i* > ieff cannot be detected (Gould 1995). The luminosity at a certain frequency u 
is Lj, = Att'^RIB^, where B^, = {2hv'^ / c") /[ex.-g{hv / kT) - 1] is the Planck function. In the 
Rayleigh-Jeans limit, e.g., K band A ~ 2.2 /xm, the Planck function is approximated by 
Bj, {2iy^ /c^)kT, and thus the luminosity in the K band is related to the stellar radius by 
R^ oc (L/f/T)^/^. Here the frequency is regarded constant over a narrow band. Since the 
temperature of relevant source star varies within a narrow range, 3 — 6 x 10^ K, relative 
to the luminosity, 1 — 10^ Lq, the dominant factor in determining stellar radius is the 
luminosity. Therefore, one can approximate the stellar radius by i?* = lO^-^^'^^'®"'^^'*^ Rq, 
where the solar K-band luminosity is Mk,q = 3.3 mag and iT-band LF are provided from 
the model LF discussed in S 6.2. 



6.2. Observational Strategy 



I begin by analyzing the event rate in a relatively modest observing program. The 
experiment is carried out using a 2K x 2K CCD camera with a pixel scale of 0".25 on 
aim telescope. If the pixel scale is too big, the geometric alignment of the images 
would introduce considerable additional noise (pixelization error; Gould 1996). I assume 
that the camera can detect 12 photons s^^ for an / = 20 star. The observation is 
carried out 4 hrs/night during the M31 season (~ 1/3 yr). Bad weather will influence 
the event rate, e.g., by degrading the S/N and by missing portions of the light curve. 
However, the weather pattern is unpredictable and differs dramatically from site to 
site. Hence, I set the observation time less than what might be being carried out 
(~ 5 — 6 hrs/night) to compensate the loss of T by the bad weather instead of modeling 
the weather pattern. The total exposure time is subdivided into smaller integration 
times, ~ 20 min, so that the images are not saturated. The angular size of the PSF is 
^psf ~ 7r0see/ln4 = 2.27 arcsec^, assuming typical seeing of ^see ~ The fractional 
seeing difference is kept to (50scc/^scc) ~ 0.05 by using methods discussed in § 7. The 
surface brightness distribution is based on the density models of the M31 bulge and disk 
and then normalized to the values of Walterbos et al. (1987). The LF is modeled primarily 
by adopting that of Galactic bulge stars shown by dotted line in the second panel of Figure 
3 (J. Frogel, private communication) 
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However, the LF in the faint end, F < Fc, is incomplete. Here Fc corresponds to 
/ = 28.2 where the LF begins to dechnc. If the LF is adopted without corrections, the 
fraction of bright stars would be overestimated. The net effect is then the systematic 
overestimation of event rate because events are more likely to be detected for brighter 
stars. Thus, I model the faint part of the LF {F < Fc) by adopting the LF of stars in the 
solar neighborhood determined by Widen, Jahrciss, & Kriiger (1983) for 3.6 < M/ < 8 
and by Gould, Bahcah, & Flynn (1996) for 8 < Mj < 14. The corrected model LF yields 
the fluctuation magnitude of Y^i 4'{Po,i)Po i/ 4'{Po,i)Po,i = 23.48 mag, which matches 
well with Tonry's (1991) determination of fiij = 24.3 before the extinction correction. The 
fraction of light from the stars in the corrected part of the LF (M/ > 3.6) is 14.7%. The 
model LF is shown with a solid line in the second panel of Figure 3. 

Combining all this information, the values of /3max are computed for a given stellar 
flux F by solving equation (2.4) and taking the additional noise term r) into consideration. 
The determined /3max as a function of the stellar flux is shown in the upper panel of Figure 
3 for the position {x',y') ~ (1.0,1.0) kpc at which ~ 16 days, E ~ 19 mag arcsec"^, 
and r ~ 4 X 10"'^ for the various (S'/A^)iiiin noted in the figure. The LF is normalized by 
the surface brightness of the position. Note that as a star becomes fainter, it can only 
be detected at higher amplification and that the detection limit decreases as (5/A^)min 
increases. Now one is ready to compute event rate for any adopted {S/N)iain with the 
results of the r and te distributions determined in § 3 and § 4. 



6.3. Event Rates 



The term "event rate" is somewhat subjective because its value depends on the 
required precision and strategy of the experiment. To measure the time scale with better 
precision, higher S/N is required. The relation between the precision of measurement 
and the required S/N is quantified by a factor A = {S/N)6te/te, where 6te/te is the 
fractional accuracy of the time scale measurement. Its value is computed as a function 
of (3 and shown in Figure 1 of Gould (1996). To determine {S/N)^[^, I use the mean 
value, (A) ~ 16, since A is a slowly varying function of (3 except in the region of very close 
encounter event (i.e., /3 < 0.01); 10 < A < 20 in the range 0.1 < P < 0.8. The required 
minimum signal-to-noise ratio is then, 

{S/NUin ~ (6.3.1) 

First I estimate the event detection rate based on the precision that can distinguish 
lensing events from flux variation due to statistical noise. In principle it is possible to 
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select events up to the Poisson limit of Ta fluctuations, i.e., (iS/A^jmin = 7 (Gould 1996). 
To be conservative, I use (S'/A^)inin = 10. The resultant event rate map dT{x' ,y')/dQ, 
dQ. = 10^ arcsec^, at this (S'/A^)min is constructed and shown in Figure 4 as a contour map. 
The contour levels are drawn in units of events yr~^. Then the total event rate detected 
over the entire area of the image, ^ccd ~ 70 arcmin^, is computed by 



resulting in 183 events yr . 

Now I stiffen the qualification of events by counting only those that can provide time 
scales with an acceptable uncertainty, 5te/te < 0.2. The required signal-to-noise ratio is 
then {S/N)iam = 80 determined from equation (6.3.1). The result is F ~ 15 events yr"-*^. 
The total event rate for other values of {S/N)inin are shown in Figure 5 for individual lens 
population events. 



There are several ways of improving the event rate. First, one can increase F by 
reducing the seeing difference, 69scc, which causes additional noise r/psf. The easiest way 
to reduce dOsee/dsee is classifying images (there will be > 1000 yr~^ of them) with similar 
seeing, and then subtracting the current image from reference images of similar seeing. 
The other method, which is being developed by Tomaney Sz Crotts (1996), is reducing the 
seeing differences by convolving the current image so that the resultant PSF matches the 
reference frame. Second, the event rate would increase by carrying out the experiment at 
more than one place. By observing more often, one could detect events with shorter t^s 
than could not be detected from a single site. The other effect of multi-observation on F is 
the increase of the time averaged photon detection rate, a. Finally, a bigger telescope will 
detect more events. The event rates with the additionally adopted observational strategies 
(single or combined) are shown in Table 1. The values in the second and third columns 
represent F when {S/N)^[n = 10 and 80, respectively. The values in the parentheses 
represent the improved factors in F compared to the value obtained under the observational 
condition described in § 6.2. 

If all of these strategies are adopted, the event rate would increase by up to a factor 
~ 6.5, so one could obtain events with time scale data at a rate ~ 62 yr~^. Therefore, a 
few years of observation would provide one with a data set large enough for the analysis of 
M31 stars and MACHOs. From this analysis, one can obtain important information such 




(6.3.2) 



7. Improved Event Rate 
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as the LF of M31 bulge stars (Gould 1996), the optical depth distribution of MACHOs, 
and above all the mass spectrum of M31 MACHOs (Han & Gould 1996a). 



8. Uncertainties 

The determined event rate is subject to many sources of uncertainties. These include 
insufficiently detailed density and velocity distributions and mass spectrum of M31 bulge 
MACHOs. Additional uncertainties from the crude models of extinction and the LF also 
propagate to the determination of the event rate. 

Among all the uncertainties, however, the most important comes from the stellar mass 
density and the fraction of dark matter in various forms, which can only be determined 
by the comparison of the theoretical predictions with the observational results. The other 
uncertainties will not seriously affect the estimates given here. For example, the velocity 
field is spectroscopically well constrained. Although the model LF at the faint end is very 
uncertain, stars in this regime make little contribution to the event rate because they are 
generally too faint to be detected. Thus, the detailed structure of the faint end of the LF is 
not important. Of course, my determination of the amount of extinction, /y, is not precise 
enough for other type of studies such as the construction of an extinction map which 
would require very detailed observation. However, in the optically thin limit, the event 
rate depends only on the mean extinction (and even for a totally opaque disk only falls by 
a factor ~ 1/2), so the error induced by uncertain extinction is unlikely to exceed 10%. 

In addition, it has been demonstrated by Tomaney & Crotts (1996) that systematic 
errors can be brought under very good control. For example, they measure the flux variation 
of a variable found in the central M31 bulge (surface brightness R ~ 18.2 mag arcsec"^) 
with a la flux error equivalent to i? = 24.6, within a factor ~ 2 of the photon limit and 8 
mag below the flux in a seeing disk. 

I would like to thank A. Gould and J. Frogel for making very helpful comments 
and kindly providing the Galactic bulge luminosity function and to USNO for providing 
telescope time. This work was supported by a grant AST 94-20746 from the NSF. 
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FIGURE CAPTION 

Figure 1: The resulting optical depth distribution caused by M31 population (M31 
disk+bulge self-lensing plus M31 halo) is shown as a function of projected position (x',y'). 
The total optical depth both by M31 and Galactic halo lenses are obtained by adding 
TMW,haio = 4.4 X 10"^ to the values marked on the contour. The contour levels are marked 
on the map in units of 10~^. The lower half of the figure contains the near side of the M31 
disk. 

Figure 2: Median time scale distribution as a function of (x', y') for M31 self-lensing 

events. The contour levels are drawn in units of days. Other population events, Galactic 
halo and M31 halo events, the time scale distributions are similar to that of M31 self-lensing 
events. 

Figure 3: Values of /3max as a function of luminosity of M31 bulge stars (upper panel) and 
the adopted LF of M31 bulge stars (lower panel). For the computation of /3max5 values of S 
and Lo at {x',y') = (1.0, 1.0) kpc are used. The basic LF function (dotted line) is modeled 
based on the Galactic bulge LF (J. Frogel, private communication). The faint end of the 
LF (solid line) is modeled by combining the LFs of Gould et al. (1996) and Wielen et al. 
(1983). 

Figure 4: The event rate per angular area, dO, = 10^ arcsec^, as a function of position 
{x', y') for the threshold signal-to-noise ratio of {S/N)-aiin = 10. The contour levels are 
drawn in units of events yr~-^. 

Figure 5: The total event rate, F = {dV /d^)Qccd: for individual population events as a 
function of {S/N)inm- 
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improved 
strategies 


{S/NUn = 10 


r 

(5/iV)„,in = 80 


Se^ee/Osee = 0.05 

2-m telescope 


366 (2.0) 


43.7 (2.8) 


Se,ee/0,ee = 0.05 

three-site observation 


293 (1.6) 


31.2 (2.0) 


Se^ee/Osee = 0.01 


201 (1.1) 


18.7 (1.2) 


69,ee/0,ee = 0.01 

three-site observation 


384 (2.1) 


43.7 (2.8) 


^^^scc/^scc = 0.01 

2-m telescope 


403 (2.2) 


48.3 (3.1) 


Se^ee/Osee = 0.01 


732 (4.0) 


101.4 (6.5) 



three-site observation 
2-m telescope 



Table 1: The event rates with improved observational strategies in units of events yr~ . 
The values in the parentheses represent the improved factors compared to that under 
observational condition described in § 6.2. The first column shows the additionally 
adopted strategies, and the values in the second and third columns are those of F when 
(5/A^)min = 10 and 80, respectively. 
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